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Bovine coronavirus (BCoV) is found worldwide and causes respiratory infections and diarrhoea in calves 
and adult cattle. In order to investigate the molecular epidemiology of BCoV, 27 reverse transcription 
polymerase chain reaction (RT-PCR) positive samples from 25 cattle herds in different parts of Sweden 
were analysed. A 1038-nucleotide fragment was PCR amplified and directly sequenced. 

The analysed BCoV strains showed a high sequence identity, regardless of whether they were obtained 
from outbreaks of respiratory disease or diarrhoea or from calves or adult cattle. Circulation of an iden¬ 
tical BCoV strain during a 4-month period was demonstrated in calves in one dairy herd. In a regional 
epizootic of winter dysentery in Northern Sweden, highly similar BCoV strains were detected. In the 
Southern and Central regions, several genotypes of BCoV circulated contemporaneously, indicating that 
in these regions, which had a higher density of cattle than the Northern regions, more extensive trans¬ 
mission of the virus was occurring. Identical BCoV sequences supported the epidemiological data that 
inter-herd contact through purchased calves was important. Swedish BCoV strains unexpectedly showed 
a high homology with recently detected Italian strains. This study shows that molecular analysis of the 
spike (S) glycoprotein gene of BCoV can be a useful tool to support or rule out suspected transmission 
routes. 

© 2012 Elsevier Ltd. All rights reserved. 


Introduction 

Bovine coronavirus (BCoV) belongs to the family Coronaviridae, 
order Nidovirales and has a single-stranded, positive-sense RNA 
genome 31 kb in length, including 13 open reading frames (Clark, 
1993; Saif, 2004). Variations in host range and tissue tropism of 
coronaviruses are attributed to the spike (S) glycoprotein (Galla¬ 
gher and Buchmeier, 2001 ). This protein is cleaved by an intracel¬ 
lular protease into two functional domains (Gallagher and 
Buchmeier, 2001). The SI subunit is peripheral and responsible 
for virus binding to host-cell receptors (Kubo et al., 1994), induc¬ 
tion of neutralizing antibodies (Yoo and Deregt, 2001), and hae- 
magglutination activity (Schultze et al„ 1991). The SI sequence is 
variable, and mutations in this region have been associated with 
altered antigenicity and virus pathogenicity (Fazakerley et ah, 
1992; Ballesteros et ah, 1997). The sequence of the S2 subunit is 
more conserved and this subunit is responsible for cell membrane 
fusion activity (Yoo et ah, 1991). The SI gene has been exploited as 


* Corresponding author. Tel.: +46 18 671782. 

E-mail address: madelein.traven@slu.se (M. Traven). 

1090-0233/$ - see front matter © 2012 Elsevier Ltd. All rights reserved. 
http://dx.doi.Org/10.1016/j.tvjl.2011.12.015 


a target to study the molecular epidemiology of BCoV infection 
(Hasoksuz et ah, 2002a; Jeong et ah, 2005; Liu et ah, 2006). 

BCoV infection is found worldwide and causes calf diarrhoea 
(CD), winter dysentery (WD) in adult cattle, and respiratory infec¬ 
tions in cattle of all ages (Clark, 1993; Paton et ah, 1998; Saif, 
2004). Infection has a high morbidity but a low mortality (Alenius 
et ah, 1991; Saif, 2004). Outbreaks typically occur in autumn and 
winter, probably in association with housing, but severe outbreaks 
in adult cattle have been reported in summer (Park et ah, 2006; 
Decaro et ah, 2008a). Economic losses can be heavy due to a 
marked reduction in milk yield (Traven et ah, 2001; Saif, 2010). Calf 
diarrhoea and respiratory disease also cause significant losses 
(Svensson et ah, 2003; Torsein et ah, 2011), but the role of BCoV 
in these diseases has not been fully evaluated. 

Several serological studies have shown that BCoV infections are 
highly prevalent in Swedish cattle herds, causing disease in both 
calves and adults (Traven et ah, 1999; Hagglund et ah, 2006; Bido- 
khti et ah, 2009). In a previous study we showed that all cattle sam¬ 
pled during a herd outbreak of BCoV were carrying identical BCoV 
sequences (Liu et ah, 2006). The objective of the present study was 
to investigate the molecular epidemiology of BCoV by analysing 
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the polymorphic region of the S gene from herds with outbreaks of 
diarrhoea or respiratory disease and to determine if this method 
could be used to support or rule out suspected transmission routes 
between herds. 

Materials and methods 

Samples 

Faecal and nasal swab samples were collected from housed cattle in 23 dairy 
herds and two feedlot herds (Table 1) in different parts of Sweden over a 4 year per¬ 
iod. Samples were collected October 2005 to May 2006, October 2006 to May 2007, 
October 2007 to May 2008, and October 2008 to July 2009. Sampled cattle in all 
herds were showing clinical signs of BCoV infection and the presenting symptoms 
are summarised in Table 1. 

The old Swedish BCoV strain, Cl-9202, which was detected in 1992 from a WD 
outbreak in Uppland and previously used to experimentally reproduce disease 
(Traven et al., 2001; Niskanen et al., 2002), was also sequenced in this study for 
phylogenetic comparison. Data on BCoV vaccination status was not available for 
the herds sampled in this study, but such vaccine is rarely used in Swedish cattle 
herds. 

RNA extraction, synthesis of cDNA and PCR 

RNA extraction with TRIzol LS reagent (Invitrogen) and cDNA synthesis with 
random priming was performed as described by Liu et al. (2006). Polymerase chain 
reaction (PCR) was performed as described in that paper with minor modifications. 
Briefly, 5 pL of cDNA sample was added to a tube containing 45 pL of PCR mixture, 
comprising 5 pL of 10 x PCR buffer, 1 pL of 10 mM dNTP mix, 5 pL of 1 mg/mL bo¬ 
vine serum albumin, 1.5 pL of each primer (10 pM), 5 pL of 25 mM MgCl 2 , 1 U of 
Taq DNA polymerase (AmpliTaq; Perkin-Elmer) and 24 pL of water. Two primers, 
S1AF and S1CR (Hasoksuz et al., 2002b) were used to amplify a 1550-nucleotide 
(nt) fragment using Biometra T3000 Thermocycler (LABREPCO). The thermocycling 
profile included 35 cycles of denaturation at 94 °C for 45 s, annealing at 50 °C for 
60 s, and extension at 72 °C for 200 s, and a final extension at 72 °C for 7 min. 

Nested and semi-nested PCR (N- and SN-PCR) assays were developed to in¬ 
crease the sensitivity of the PCR method for clinical detection by amplifying 5 pL 
of the PCR products under the same reaction conditions described above, but with 
a 120 s extension time, to amplify the target fragments using published primers 
(Hasoksuz et al., 2002b); S1AF/S1AR primer pair for the initial 654-nt fragment 
(fragment A) as SN-PCR, S1BF/S1BR primer pair for the second 490-nt fragment 
(fragment B) as N-PCR. 


DNA sequencing and phylogenetic analysis 

The PCR products were purified and subsequently sequenced using the same 
primers as for PCR and an ABI PRISM BigDye Terminator v3.1 Cycle Sequencing 
Kit (Applied Biosystems) as previously described (Liu et al., 2006). Sequencing 
was performed at least twice in each direction in an ABI 3100 genetic analyser (Ap¬ 
plied Biosystems). The sequencing reads for each sample were assembled to a final 
1038-nt fragment of the SI subunit spanning nucleotides 23641-24679 (amino 
acid residues 1-346 of the S glycoprotein) of BCoV strain Mebus (GenBank acces¬ 
sion number U00735). Multiple sequence alignment was performed with software 
BioEdit version 7.0.9.0 (Hall, 1999) with an engine based on the ClustalW algorithm. 

Maximum likelihood analysis was performed with software PHYML (Guindon 
and Gascuel, 2003), as previously described (Liu et al., 2009). The General Time 
Revisable (GTR) with substitution rate heterogeneity and proportion of invariable 
sites ‘GTR +1 + G’ was selected for inferring BCoV phylogeny, and 1000 non-para- 
metric bootstrapping analyses were performed to evaluate the robustness of the 
tree (Guindon and Gascuel, 2003). The dataset included the reference strains of 
BCoV: Mebus (GenBank accession number U00735); Quebec (AF220295); Ger- 
many/V270 (EF193075); alleged respiratory bovine coronavirus strains: LSU2 
(AF058943) and OK-3 (AF058944); alleged enteric bovine coronavirus strains: FI 5 
(D00731) and LY138 (AF058942). The analysis also included recently published se¬ 
quences of three Italian strains 179/07-11 (EU019216), 438/06-TN (EU814647), and 
339/06 (EF445634); two Korean strains: KWD19 (DQ389660) and KCD10 
(DQ389641), and one Japanese strain: Kakegawa (AB354579). One Swedish strain 
RclN-99 (DQ121649) and two Danish strains AclF-03 (DQ121619) and HclF-01 
(DQ121628) from the previous study (Liu et al., 2006) were included in the analysis, 
where 414 nt were treated as missing data. The analysis also included Human cor- 
onaviruses: Hu-OC43 (Z32769) and Hu-4408 (L07748); Giraffe coronavirus 
(EF424624) and Alpaca coronavirus (DQ915164). 


Results 

Genetic variation of the S gene 

Comparative analysis of the 1038 nt sequence showed that all 
the 26 Swedish strains obtained in 2005-2009 shared a high iden¬ 
tity both at the nt level (>97.7%) and at the deduced amino acid 
(aa) level (>96.7%). Compared with the reference strain Mebus, 
29-36 nt substitutions were identified, resulting in 16-20 aa 
changes that were clustered around residues 11-40, 88-115, 
143-188, and 241-257. 


Table 1 

Location of herds positive for BCoV, sampling time, animals sampled, symptoms and type of sample. 


Herd 

Type 

Year/month 

Age (symptoms b ) 

Sample type 

County 

AC1 

Daily 

2006/11 

Adult (WD) 

Faecal 

Vasterbotten 

AC2 

Daily 

2007/11 

Adult (WD) 

Faecal 

Vasterbotten 

Cl 

Daily 

1992/02 

Adult (WD) 

Faecal 

Uppland 

C2 

Dairy 

2007/01 

Adult (WD) 

Faecal 

Uppland 

C3 

Dairy 

2007/11 

Adult (WD) 

Faecal 

Uppland 

C4 

Dairy 

2007/12 

Adult (WD) 

Faecal 

Uppland 

C5 

Dairy 

2008/01 

Adult (WD) 

Faecal 

Uppland 

C6 

Feedlot 

2009/03 

Calf (RD) 

Nasal 

Uppland 

11 

Daily 

2007/02 

Adult (WD) 

Faecal 

Gotland 

12 

Daily 

2007/02 

Adult (WD) 

Faecal 

Gotland 

13 

Dairy 

2007/03 

Adult (WD) 

Faecal 

Gotland 

14 

Dairy 

2008/10 

Adult (WD) 

Faecal 

Gotland 

Ml 

Dairy 

2006/02 

Calf (CD) 

Faecal 

Skane 

M2 a 

Dairy 

2006/05 

Calf (CD) 

Faecal 

Skane 

N1 

Daily 

2005/11 

Calf (CD) 

Faecal 

Halland 

N2 a 

Daily 

2006/03 

Calf (CD) 

Faecal 

Halland 

N3 

Daily 

2006/05 

Calf (CD) 

Faecal 

Halland 

N4 

Dairy 

2009/03 

Adult (WD) 

Faecal 

Halland 

PI 

Dairy 

2009/02 

Adult (RD) 

Nasal 

Vastergotland 

U1 

Feedlot 

2009/07 

Calf (RD) 

Nasal 

Vastmanland 

Y1 

Dairy 

2008/01 

Adult (WD) 

Faecal 

Angermanland 

Z1 

Daily 

2007/11 

Adult (WD) 

Faecal 

Jamtland 

Z2 

Daily 

2007/11 

Adult (WD) 

Faecal 

Jamtland 

Z3 

Daily 

2007/11 

Adult (WD) 

Faecal 

Jamtland 

Z4 

Dairy 

2007/11 

Adult (WD) 

Faecal 

Jamtland 

Z5 

Dairy 

2007/11 

Adult (WD) 

Faecal 

Jamtland 

Z6 

Dairy 

2007/12 

Adult (WD) 

Faecal 

Jamtland 


a Herds Ml and N1 visited on separate occasions, designated M2 and N2. 
b WD, winter dysentery; CD, calf diarrhoea; RD, respiratory disease. 
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More genetic variation was observed within Southern strains 
(23 nt changes) than within Northern and Central strains (9 nt 
changes). Sequences from five herds (Z1-Z5) in the Northern re¬ 
gion (Jamtland) sampled during the third year were identical. 
These sequences were also identical to the sequence Ul-0907 ob¬ 
tained from a herd in the Central region (Vastmanland) that had an 
outbreak of respiratory disease in July 2009. The two sequences 
Nl-0511 and N2-0603 obtained from the same herd with a 4- 
month interval showed 100% identity, while 6 nt changes were 
found between sequences Ml-0602 and M2-0605 obtained with 
a 3-month interval from another herd. 

The BCoV strains obtained from respiratory and enteric disease 
did not show any consistent differences in the sequenced SI re¬ 
gion. Strain Cl-9202 shared a 97.0-98.2% identity at the nt level, 
96.2-97.4% at the deduced aa level, with the Swedish strains de¬ 
tected in 2005-2009. GenBank accession numbers for the BCoV se¬ 
quences obtained in this study are JN795141-JN795167. 


Phylogenetic analysis 

The Swedish BCoV strains obtained for this study fell into four 
clusters in the consensus tree based on maximum likelihood anal¬ 
ysis of the 1038-nt fragment of the S gene (Fig. 1). Cluster I in¬ 
cluded 20 strains from all over Sweden: three from Halland and 
one from Vastergotland in the South; 2 from Gotland, an island 
in the east; 5 from Uppland and 1 from Vastmanland in the Central 
region; 6 from Jamtland and one each from Vasterbotten and 
Angermanland in the North. Cluster II contained 4 Swedish strains: 
AC1-0611 from the Northern region Vasterbotten and Ml-0602, 
M2-0605 and N4-0903 from the Southern regions Skane and Hal¬ 
land, and a Danish strain AcF-03 (Liu et al„ 2006). This cluster also 
included the Italian strain Bubalus/ItaIy/179/07, a BCoV strain de¬ 
tected in buffalo calves with severe enteritis (Decaro et al„ 2008b). 
Cluster 111 contained the two Swedish strains 12-0702 and 14-0801 
from Gotland and the strains Italy/438/06 and Italy/339/06 that 



Z30711 

Z2-0711 

Z1-0711 

Ul-0907 

Z5-0711 

Z4-0711 

C5-0801 
C4-0712 
Z6-0712 
Y1-0801 
AC2-0711 
L-C34J711 
N3-0605 
PI-0902 

C6-0903 

N2-0603 

Nl-0511 

134703 

C2-O701 

11-0702 

M2-0605 

Ml-0602 
AclF-03 

74 | N4-0903 

AC1-0611 

Bubalus/ltaty/179/07 


Alpaca 


■ Hu-OC43 


P.0050 


N 



A 


B 


Fig. 1. A majority-rule consensus tree (A) of maximum likelihood analysis of the Swedish BCoV strains sampled across the country (B) during four in-stable seasons. Four 
clusters are labelled I, II, III, and IV, respectively. Bootstrapping values (italic) are placed next to the major clades. Bar indicates changes per site. The tree is presented as rooted 
with LY-138, FI5 and Hu-OC43. Herd Cl is not included on the map. 
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were isolated from an outbreak of enteric and respiratory disease 
in dairy cattle (Decaro et al., 2008a). Cluster IV contained Cl- 
9202, the oldest Swedish strain from Uppland, together with the 
Swedish strain RclN-99 and the Danish strain HcF-01 that were 
collected in 1999 and 2001, respectively. The grouping of these 
three sequences was supported by a bootstrapping value of 62%. 


Discussion 

Twenty-seven samples from 25 herds located in different parts 
of Sweden were positive for BCoV and comparative sequence anal¬ 
ysis showed a high identity between the strains. There were no 
particular aa residues of the initial 346 aa of the SI subunit that 
could differentiate strains obtained from cattle with respiratory 
disease from those obtained from calves or adult cattle with diar¬ 
rhoea. The results suggest that the same strains of BCoV cause nat¬ 
ural outbreaks of respiratory disease, CD and WD, which is in 
agreement with previous studies (Traven et al„ 2001; Niskanen 
et al., 2002; Saif, 2010). 

The occurrence, in this study, of a BCoV outbreak in summer 
confirms previous reports from South Korea and Italy (Park et al., 
2006; Decaro et al., 2008a) that BCoV circulates between cattle 
herds all the year round. Furthermore, the identical 1038-nt se¬ 
quences of strain Ul-0907, obtained from an outbreak of respira¬ 
tory disease in calves, and strains Z1-0711-Z5-0711, obtained 
from WD outbreaks in Jamtland, indicate that one BCoV strain 
may circulate in the cattle population for more than a year. Herd 
U1 is a feedlot which regularly purchases calves from Jamtland. 
Analysis of the whole S gene (4092 nt) of strains Z2-0711 and 
Ul-0907 confirmed an identity of 99.6%, suggesting that sequence 
analysis of the S gene can be used to support or rule out suspected 
transmission routes of BCoV between regions. 

The finding of identical 1038-nt sequences of two strains ob¬ 
tained 4 months apart from young calves in a large dairy herd 
(Nl, >200 cows) indicates, for the first time, that a BCoV strain 
may circulate for an extended period of time within large herds. 
On the other hand, strains from another large herd (Ml) obtained 
3 months apart had several mutations, suggesting introduction of a 
new strain into that herd. 

All strains obtained from Northern and Central Sweden during 
the third year were highly similar (nt identity >99.7%). These 
strains were obtained from a regional epizootic of WD, i.e. a series 
of outbreaks close in time and space. These regions were shown to 
have a low herd immunity to BCoV; 70% of the herds had BCoV 
seronegative primiparous cows in the Northern regions compared 
to 10% in the Southern parts in 2007 (Ohlson et al., 2010; A. Ohlson 
et al., unpublished data). This may be related to a lower density of 
herds in Northern regions, less transportation of animals, and/or a 
difference in biosecurity between the regions (Ohlson et al., 2010). 
The low herd immunity suggests that BCoV is not continuously cir¬ 
culating in the Northern parts of Sweden, and that virus causing re¬ 
gional epizootics is introduced from other regions. Such a 
transmission is also suggested by the finding of strain N3-0605 
in Southern Sweden prior to and highly identical (nt identity 
>99.8%) to the strain causing the Jamtland outbreaks of WD (Zl- 
Z5). 

The molecular pattern of clusters II, III and IV suggests that the 
BCoV strains circulating in Swedish herds are closely related to 
strains present in other parts of Europe. While it is easy to under¬ 
stand the clustering of strains from the Central and Southern parts 
of Sweden with strains from Denmark, which borders with Swe¬ 
den, grouping of Swedish strains with Italian strains was surpris¬ 
ing. As there were no records of animals purchased from Italy, it 
is unlikely that the virus was introduced to Sweden by transporta¬ 
tion of animals directly from Italy. There have, however, been ani¬ 


mals exported from Sweden to Central and Southern Europe in 
recent years. Samples sequenced from other European countries 
are needed in order to trace the origin and/or spread of Swedish 
strains and molecular analysis of future BCoV outbreaks from dif¬ 
ferent parts of the world will increase the understanding of BCoV 
epidemiology, e.g. direct and indirect transmission routes between 
herds, regions or even countries. 

Conclusions 

In this study we showed that molecular analysis of the S gene of 
BCoV was a useful tool for studying the epidemiology of this virus. 
The analysed strains showed a high sequence identity, regardless 
of whether they were obtained from outbreaks of respiratory dis¬ 
ease or diarrhoea or from calves or adults. Circulation of a BCoV 
strain during an extended period of time was demonstrated in 
calves in a large dairy herd. In Northern Sweden, a regional epizo¬ 
otic of highly similar BCoV strains was detected, while in the 
Southern and Central regions more extensive transmission of the 
virus was indicated. Finally, identical BCoV sequences supported 
epidemiological data on inter-herd contact through purchased 
calves. 
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